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ABSTRACT 

Axion-likc particles (ALPs) belong to a class of new pseudoscalar particles that generically couple to 
photons, opening the possibility of oscillations from photons into ALPs in an external magnetic field. 
These oscillations are expected to imprint irregularities in a limited energy range of the spectrum 
of astrophysical sources, witnessing the turbulence of their magnetic fields. In this study, Chandra 
observations of Hydra galaxy cluster are used to constrain the value of the coupling of ALPs to photons. 
The case of the conversion of X-ray photons from the central source Hydra A in the magnetic field of the 
cluster is considered. The magnetic field strength and structure is well determined obscrvationally, 
which adds to the robustness of the analysis. The absence of anomalous irregularities in the X- 
ray spectrum of Hydra A conservatively provides the most competitive constraints on the coupling 
constant for ALP masses below 7 x 10~^^ eV at the level of g^a < 7.3 x 10^"'^^ GcV^^. Because of the 
specific phenomenology involved these constraints actually hold more generally for very light pseudo 
Nambu-Goldstone bosons. 

Subject headings: Elementary particles — Galaxies: clusters: individuahHydra A — X-rays: individ- 
ual: Hydra A 



1. INTRODUCTION 

Many extensions of the standard model of particle 
physics predict the existence of new light bosons. One 
of the most studied example is the axion, a pseudoscalar 
particle associated to the spoi itaneous breakiiig of th e 
U(l) Peccei-Quinn symmetry (jPeccei &: QuinnI [l977b[ ). 
The symmetry was introduced as a solution to the 
strong CP problem, with the spontaneous breaking of 
the symmetry dynamica lly tuning the CP angle to zero 
(jPeccei fc QuinnI ri977aD . In this case, the mass of the 
axion is predicted to scale with the inverse symmetry 
breaking scale /. Pseudoscalar particles generically cou- 
ple with the electromagnetic field via a two photon vertex 
of the form: 

Cja = --g^aF^^F^^a = g^aE ■ Ba , (1) 

where F is the electromagnetic tensor, F its dual, E and 
B the electric and magnetic fields, a the axion field and 
g-ya the axion-photon coupling strength. For standard 
axions, g-ya also scales with 1// so that the coupling 
strength and the mass of the hypothetical particle are 
proportional. However pseudoscalar particles with unre- 
lated mass and coupling strengt h can be e xpected from 
more general U(l) symmetries fKiin[ JT987l ) or in extra- 
dimension gauge theories ( JTurok 1996| ). Those have the 
same phenomenology than standard axions and are called 
axion-like particles (ALPs). The coupling of ALPs to 
two photons implies the possibility for oscillatio ns from 
photons into ALPs in an external magnetic field (jSikivid 
[1983: Raffch & Stodolsky 1988). This is used by labora- 
tory experiments to search for ALPs produced in the Sun 
(|Arik et al.l 120091 ) or by conversion of a laser beam in a 
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strong magnetic field (|Ehret et al.ll2009() , to put stringent 
constraints on the coupling strength. 

Astrophysical environments offer the possibility of 
strong magnetic fields on long baselines and are thus 
promising targets in the search for ALPs. In contrast 
with laboratory experiments, magnetic fields in astro- 
physical environments are usually not coherent. In the 
case of galaxy cluster magnetic fields for instance, the 
turbulence is well described by a Kolmogorov power 
spectrum on scales from tens of par secs up to a few 
kilopa rsecs. It has been shown in iWouters fc BrunI 
(|2012l ) that when 7-rays mix with ALPs in a turbu- 
lent magnetic fields, the turbulence of the field trans- 
lates into an irregular behavior of the 7-ray energy spec- 
trum. These irregularities are expected in a limited 
energy range, for energies around the energy thresh- 
old of the mixing. The search for irregularities in the 
TeV-energy spectrum of the bright blazar PKS 2155- 
304 measured by H.E.S.S. gives stringent constraints 
in a limited range for the ALP mass around 1 0~^eV 
(|Wouters fc Brun for the H.E.S.S. collaboration! |2013| ). 
In the present study, X-ray data are used and ALP 
masses below the effective mass of the photon in the 
medium are considered. For typical values of the elec- 
tron density in galaxy clusters, the effective photon mass 
is about 10~^^ eV. Spectral irregularities associated to 
this low mass are expected in the X-ray band for typical 
galaxy cluster magnetic fields and values of the coupling 
g^a not excluded by current constraints. One possible 
source of uncertainty and systematic error in the deter- 
mination of constraints on ALPs in astrophysics is the 
lack of knowledge about the magnetic fields that serve 
as a target to trigger the photon/ ALP conversions. Here 
we choose to use one of the best characterized large scale 
magnetic field in both strength and structure that is the 
magnetic field of the well studied Hydra cluster. In the 
following, irregularities are searched in the spectrum of 
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the bright X-ray source Hydra A that lies at the center 
of the galaxy cluster Hydra. The characteristics of the 
electron density and those of the magnetic field (profiles, 
scales, turbulence) in the cluster derived by Faraday ro- 
tation measure have been extensively discussed in the lit- 
erature and are used in the analvsis (fVoat fc Enfilin 20031 
[20051: iKuchar fc Ei"Imml[20Tlt iLaing et al.ll2008D . The ar- 
ticle is organized as follows. First, the phenomenology 
of the photon/ ALP oscillations in the medium is briefly 
recalled. In a second step, the modeling of the magnetic 
field and electron density in the cluster is presented, fol- 
lowing the most recent measurement by Faraday rota- 
tion. The results of the analysis of Chandra data on Hy- 
dra cluster are then exposed and the derived constraint 
is finally discussed. 

2. PHENOMENOLOGY OF THE PHOTON/ALP SYSTEM 

The photon/AL P sys tem is described as in 
iRaffelt fc Stodolskvl ()1988() by a wave-function with 
three states, two for the photon corresponding to 
the polarization states and one for the axion. The 
Lagrangian of Eq. [T] induces the mixing between the 
three states in an external magnetic field. A transfer 
matrix is used to compute the probability of observing 
a photon after propagation in a magnetic domain of 
size s, starting from an unpolarized beam of photons. 
This probability is called the survival probability. For 
propagation in one domain with a coherent magnetic 
field of strength projected on the polarization plane B, 
this probability writes: 

where E is the energy of the photon and Ec = \m'i — 

m'^\/2gjaB defines the critical energy above which the 
mixing is efficient, m^ = Anane/mi, is the effective mass 
of a photon propagating in a plasma with electron den- 
sity rig. In galaxy clusters, typical values for the elec- 
tron density are Ue ~ 0.01 cm~'^ corresponding to an 



effective photon mass tti. 



3 X 10"" eV. For ALP 



masses ma negligible compared to to-,,, the critical en- 
ergy no longer depends on ma- The survival probability 
is energy-dependent for energies around E,,- For turbu- 
lent magnetic flelds, the global transfer matrix of the 
system has a very complex energy dependence and the 
survival probability can show strong spectral ir r egular - 
ities around Ec. As shown in iWouters fc Bruiil ()2012D . 
the exact structure of these irregularities depends on the 
realization of the turbulent magnetic field crossed by the 
beam and is therefore not predictable. However, the sta- 
tistical properties of the spectral irregularities are a pre- 
diction of the model. The following analysis addresses 
the question of the level of spectral irregularities induced 
by photon/ALP oscillations that can be accommodated 
by the data. ALPs with a low mass, lower than the 
effective photon mass, are considered. In this regime, 
m,a "C 771-), ~ 10~^^ eV, Ec does not depend anymore on 
m.a and the constraint on g^a that is deduced does not 
depend on the mass of the ALP. In this mass region, the 
coupling strength is limited to g^a ^ 10~^^GeV~^ be- 
cause of th e non-observation of a 7-ray counterpart to 
SN 1987A (jBrockwavet al.lll996D . For this upper limit 
value and typical values of magnetic fields in galaxy clus- 



ter, _B ~^ 10 fiG, Ec lies around a few tens of keV, thus 
motivating observations in X-rays. 



3. MODELING OF THE GALAXY CLUSTER 

For one to obtain constraints on ALP parameters, 
a bright X-ray point-like source embedded in a strong 
magnetic field is required. A good knowledge of the 
magnetic field and the electron density is also essen- 
tial in order to limit possible systematic bias from 
these uncontrolled parameters. One of the best candi- 
date source is Hydra A, a Fanaroff-Riley class I radio 
galaxy located at redshift z = 0.0538 that has been 
observed by the Chandra observatory. It is centered 
on a corona of hot thermal electrons emitting in X- 
rays by radiative coolin g. The electron density pro- 
file from IVogt fc EnfilinI ()2005l ) is used. It i s estimated 
with the X-ray surface brightness model of iMohr et al.l 
(jiggj ) from ROSAT PSPC data, depro.jected of t he liii e 
of sight with the method of iPfrommer fc EnfilinI ()2004[) . 
The magnetic field profile surrounding Hydra A has 
been extensively studied using Faraday rotation maps 
of the polarization of the radio emission of the AGN 
lobes due to the propagation through t he magnetized 
electron plasma (iTavlorfc PerlevI 1199 3": Vogt & En61in| 
[20031 [200a iLaing et al.l I2008t IKuchar fc Enfilin 20 Hi) . 
These studies assume that the magnetic field scales with 
the electron density as B(r) oc ndr)"^ where as is 
a free parameter to be determined. The geometry, 
described by the angle of projection of the northern 
jet on the line of sight 9, plays a crucial role in de- 
termining the strength of the magnetic field. A de- 
polarization asy mmetry is observed b etween the two 
lobes of Hydra A (jTavlor fc Perlevlll993l). most likelv re- 
lated to the Laing-Garrington effect (jGarrington et aLl 
1988: Laing 1988), beca use of the non- vanishing 9 angle. 
Tavlor fc PerlevI (J1993D found a most likely value for 



of 45°. This value has been confirmed in more re cent 
analyses (jKuchar fc EnfilinI 120 111 : iLaing et al.ll2008| ) but 
the associated uncertainty remains la rge. For instance, 
a val ue as low as 30° is still plausible (IKuchar fc EnfilinI 
120 111) . In the analysisof'Kuchar fc Enfilinj (|201lO, which 
is a refinement of lVogt fc Enfilin ( 2005 ), the strength of 
the magnetic field Bq at the center of Hydra A is found 
to be 21 fiG a 9 = 30°, compared to Bq = 36 /iG if 
9 = 45°. For cons ervativeness, the p rofile of the mag- 
netic field found in IKuchar fc EnfilinI ([20TTI ) for 9 = 30° 
is assumed in the following. The most likely value of 
the scale parameter is as — 1, which is significantly 
higher than the theoretical value as — 0.5 expected 
if the ratio of magnetic to kinetic energy is constant 
throughout the cluster. Again, the value as = 1 is 
taken in order to be conservative. Finally, the turbu- 
lence power spectrum is found to be in good agreement 
with a Kolmogorov spectrum P{k) ex k~^'^ up to scales 
as large as 10 kpc. When i t comes to spectral irreg ular- 
ities, it has been shown in IWouters fc BrunI (l2012f ) that 
photon/ALP mixing rapidly becomes irrelevant for lower 
scales for two reasons. First, the decreasing Kolmogorov 
power spectrum suppresses contributions on lower scales 
as s"^/^ (jMirizzi et al.ll2007t ). Second, the amplitude of 
photon/ALP also depends on the turbulence scale, inde- 
pendently of B, as s~^/^. For these reasons, in the fol- 
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lowing the turbulence power spectrum is modeled with a 
Kolmogorov slope between scales of 1 kpc to 10 kpc. 

4. CHANDRA DATA ANALYSIS 

Hydra A has been observed by the ACIS instru- 
m.ent onboard the sat ellite borne Chandra observatory 
(jWeisskopf et al.ll2000( ) in 1999 and 2003, for respectively 
20% and 80% of the total live-time. The ACIS is com- 
posed of two arrays of imaging CCDs sensitive to X-rays 
between 0.2 and 10 keV. It features an average angular 
resolution of 1" and an energy resolu tion of about 0.1 
keV at 1.5 keV (jGarmire et al.ll2003D . The data have 
been recalibrated using the standard procedure with the 
latest calibration 4.5.5.1 and the analysis is performed 
with CIAO tools version 4.5. Events with energy between 
0.3 and 10 keV have been retained for the analysis and 
checked for background event flares. A total live-time ex- 
posure of 238 ks is eventually available for the analysis. 
Images of th e Hydra A region observ e d by Chandra has 
been g iven in iMcNamara et al] ()2000[) ; iKirkpatrick et all 
(|2009( ) and show a diffuse emission from thermal elec- 
trons surrounding the central source. In this study, only 
the non-thermal component of the spectrum from the 
central source is of interest. To extract this spectrum, an 
aperture of 1" around the position of th e central source 
(dete rmined from infrared observations (jSkrutskie et al.l 
l2006f )) is used. Contamination from the thermal back- 
ground beneath the non-thermal component is estimated 
from an annular region between 1" and 2.5 ". The 
source spectrum is rebinned to ensure a minimum of 
30 counts per bin. The spectral analysis is carried out 
with the XSPEC package version 12.7.1 and the instrumen- 
tal response functions generated via the specextract 
procedure. The spectrum is well modeled (x^/nd.o.f. 
= 48.59/58) by a power law {dN/dE ex E'^ , T = 
1.52 ± 0.17) heavily absorbed by a hydrogen column 
A^H = 2.54 ± 0.33 X 10^2 cm'^ (jMcNamara et al.ll2000l ). 
The absorbed integrated flux of the source between 2 and 
10 keV is 3.10 ± 0.07 x IQ-^ -f.cm-'^s-^. The hydrogen 
column found with the fit is hundreds of times larger 
than the Milky Way contribution expected in the direc- 
tion of Hydra A. Evidence for an opaque system at the 
core o f Hydra A have already been reported in iTavlori 
(jl996[) . Because of the heavy absorption at low energies, 
the non-thermal component from the core of Hydra A is 
only visible above 1 keV. In the following, the spectrum 
is restricted to energies above this threshold. The spec- 
trum is shown on Fig. [1] together with the model. The 
solid line is the best-fit function, corresponding to an ab- 
sorbed power-law. The small features that appear on the 
continuous line are related to absorptions due to elements 
heavier than hydrogen. Also shown are the residuals of 
the fit normalized to the model. No significant deviations 
or irregularities that could be linked to photon/ ALP os- 
cillations are observed in the spectrum so that it is now 
used to constrain the value of f;^a- 

5. THE CONSTRAINT 

To estimate the maximum level of irregularities in- 
duced by photon/ ALP oscillations that can be accom- 
modated by the data, the spectrum is fitted with the 
model of Sec. 21 multiplied by an irregularity pattern 
corresponding to photon/ ALP conversi ons. An exam- 
ple of such pattern is given in [Wouters fc BrunI ()2012D . 
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Fig. 1. — spectrum of the non-thermal component of the point- 
Hke source. Top panel: reconstructed spectrum with model. Bot- 
tom panel: residuals normalized to the model. 



The value of g^a is a new free parameter that drives the 
amplitudes of the irregularities. When g^a goes to zero, 
the spectral model without the ALP pattern as shown in 
Fig. [T] is retrieved. Because of the turbulent nature of the 
magnetic field crossed by the beam, the exact structure of 
the ALP signal is not predictable. This procedure needs 
therefore to be repeated a large number of time, for differ- 
ent random realizations of the magnetic field. Note that 
this randomness accounts for the lack of knowledge on 
the magnetic field and not its dynamics, its configuration 
appearing static over periods of the order of thousands 
of years for the smallest considered scales. The average 
profile of X^{9"/a) for arbitrary small ALP mass over var- 
ious realizations of the magnetic field is shown on Fig. [2j 
The 95% envelope of the profile is also shown. For the 
majority of realizations, g^a ^ is preferred. However, 
the 95% lower envelope exhibits a clear minimum around 
g^ya ^ 0.4 X 10~^^GeV~"'^, corresponding to a better fit 
with the ALP pattern compared to the simple absorbed 
power-law. This effect is not an indication for a signal 
but rather an expected effect. For the corresponding re- 
alizations of the magnetic field, the spectral irregularities 
that are produced overfit the poissonian noise of the data 
that comes from the finite statistic. At the minimum, 
the irregularities in the flux have the same amplitude 
than the shot noise in the bins. One can roughly esti- 
mate the value of the coupling constant for which this 
is expected. Better fits should indeed be found when 
the scale of the induced irregularities matches the natu- 
ral statistical fiuctuations. There are 30 events per bin 
in the spectrum of Fig. [1] so that the relative level of 
natural irregularity is l/\/30 ~ 20%. For photon/ ALP 
induced irregularities, the relative (5(/)/0 fiuctuations are 
computed from the average of the surviv al probability 
over many domains ([Grossman et al.ll2002[ ): 
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where B^L = J^ B^{r)dr = 6.5 x IQ-^^ GcV^ is the in- 
tegrated magnetic field profile over the line of sight and 
s = 10 kpc is the coherence length. The irregularities are 
then comparable to the statistical fluctuations for g^a of a 
few 10~^^ GcV~^, in agreement with the minimum found 
in Fig. [2] The value of g^a at this minimum is thus pre- 
dicted as a function of the magnetic field parameters and 
the statistic in the bins. For this reason, the minimum 
does not correspond to an indication for an ALP. More- 
over, there is no guarantee on the exact spectral model 
underlying the data. There could well be for instance line 
features due to the fluorescence of iron or other heavy el- 
ements that would mimic the ALP pattern. In the end, 
this method could not be used for a discovery, but gives 
conservative constraints. For higher coupling strengths, 
the irregularities become larger than the typical range of 
noise induced by the Poisson statistic and the fit degrades 
quickly. The constraint on g^a is then set using the 95% 
lower envelope. Values of g^a with Ax^ > 4, compared 
to the minimal value of the x^ are excluded. This yields 
the constraint g^a < 7.3 x IQ-^^ GeV-^ at the 95% C.L. 
This constraint is derived for ALP models with arbitrar- 
ily small masses nia- When the ALP mass compares to 
the effective photon mass in the plasma, the spectral ir- 
regularities are no longer independent of ma ■ As seen in 
Sec.[2l this occurs when nia ^ 10~^^ GcV~^. The proce- 
dure described above is thus repeated for different values 
of rUa in order to derive the shape of the limit in this 
range of mass. The contour of exclusion obtained at the 
95% C.L. is shown on Fig. |3l There is a small range of 
mass around 8 x 10^^^ cV where the exclusion is slightly 
less constraining. This corresponds to ALP masses com- 
parable to the effective photon mass in the plasma. From 
Eq. [21 the survival probability is energy independent 



when 



In this case, no spectral irregularities 



occur and no constraint can be given. However, when 
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Fig. 2. — Median profile ofx^{9ja) for arbitrary small ALP mass 
over 1000 realizations of the magnetic field (solid line). The dashed 
lines represent the 95% envelope of the profile. The grey filled area 
gives the exelusion on g-ya that is obtained from the 95% lower 
envelope. 
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Fig. 3. — Constraint on g^a from X-ray observations of Hydra A 
cluster (hatched area). The grey filled area is the constraint from 
SN 1987 A burst l IBrockwav et al.llT993) . Exclusions arc set at the 
95% C.L. 

considering the propagation in Hydra A cluster, the elec- 
tron density is not uniform and decreases along the beam 
path. At higher ALP masses, the exclusion gets less con- 
straining as the critical energy increases to energies that 
cannot be probed with Chandra. Therefore the excl usion 
curve rises as tti^ as expected. iHorns et al.l (|2012[ ) also 
sets a constraint on very low mass ALPs, but due to the 
lack of knowledge on the magnetic field properties used 
for their analysis, the authors only constrain the g^aB 
product. The constraint obtaine d from the burst of SN 
1987 A in lBrockwav et al.l (I1996D is also shown on Fig. El 
For these small masses, direct sear ches at CAST yie ld a 
limit of the order of 10"!° GeV"! lArik et al.l (pOOJ ). It 
appears that the present constraint is therefore the most 
competitive for low mass ALPs. 

A final remark is that the current limits hold for a more 
general class of particles than ALPs. Indeed they are 
valid for generic pseudo Nambu-Goldstonc bosons, even 
in the case of a F^ type of coupling to photons instead 
of the FF term introduced in Eq. |T1 Switching to the 
first type of coupling leads to swap the polarization of 
the photon that is involved in the mixing. As magnetic 
field directions are randomly chosen inside many domains 
here, the effect of the mixing on the irregularities is the 
same. 

6. DISCUSSION AND CONCLUSION 

The search for irregularities in the X-ray spectrum of 
Hydra A measured by Chandra leads to a constraint 
at the 95% C.L. on the ALP coupling to two photons 
of g^a < 7.3 X 10~^^GeV"^ for ALP masses below 
7 X 10~^^ eV. In this study, the key point is the use of a 
measured profile of the magnetic field, determined from 
Faraday rotation maps. The related measurement of the 
magnetic field yields a determination of its strength that 
has an error, which for conservativeness the lower limit 
of is retained in the present analysis. The constraints are 
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thus considered as firm and robust. For a more likely pro- 
file of the magnetic field, the constraint would improve 
to g^a < 4.3 X 10~^^GeV~^. A refined measm-ement 
of the jet orientation angle 9 of Hydra A would help to 
characterize the magnetic field profile and possibly im- 
prove the limit. The constraint derived here is limited by 
the statistic available. Simulations of spectra that would 
be observed with more data show that with a ten time 
longer exposure, the upper limit could be improved to 



g^a < 3.9 X lO^^^GeV^"'^ for the conservative magnetic 
field profile. The conservative exclusion derived in this 
work is the most stringent constraint to date in the range 
of very low-mass ALPs, with nia < 10"-^^ eV. 



This work benefited from the support of the French 
PNHE (Programme National Hautes Energies). 



REFERENCES 



Arik, E., Aunc, S., Auticro, D., et al. 2009, JCAP, 2, 8 
Brockway, J. W., Carlson, E. D., & Raffelt, G. G. 1996, 

Phys.Lctt., B383, 439 
Ehrct, K., Frede, M., Ghazaryan, S., et al. 2009, NIM A, 612, 83 
Garmire, G. P., Bautz, M. W., Ford, P. G., Nousek, J. A., & 

Ricker, Jr., G. R. 2003, in SPIE Conference Series, Vol. 4851, 

SPIE Conference Series, ed. J. E. Truemper & H. D. 

Tananbaum, 28-44 
Garrington, S. T., Leahy, J. P., Conway, R. G., & Laing, R. A. 

1988, Nature, 331, 147 
Grossman, Y., Roy, S., & Zupan, J. 2002, Physics Letters B, 543, 

23 
Horns, D., Maccione, L., Mirizzi, A., & RoncadeUi, M. 2012, 

Phys.Rev., D85, 085021 
Kim, J. E. 1987, Phys.Rept., 150, 1 
Kirkpatrick, C. C, Gitti, M., Cavagnolo, K. W., et al. 2009, 

ApJLett., 707, L69 
Kuchar, P. & Enfilin, T. A. 2011, Astron. & Astrophys., 529, A13 
Laing, R. A. 1988, Nature, 331, 149 
Laing, R. A., Bridle, A. H., Parma, P., & Murgia, M. 2008, 

MNRAS, 391, 521 
McNamara, B. R., Wise, M., Nulsen, P. E. J., et al. 2000, 

ApJLett., 534, L135 
Mirizzi, A., Raffelt, G. G., & Serpico, P. D. 2007, Phys. Rev. D, 

76, 023001 



Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517, 627 
Peccei, R. D. & Quinn, H. R. 1977a, Phys. Rev., D16, 1791 
Peccei, R. D. & Quinn, H. R. 1977b, Phys. Rev. Lett., 38, 1440 
Pfrommer, C. & Enfilin, T. A. 2004, Astron. & Astrophys., 413, 

17 
Raffelt, G. & Stodolsky, L. 1988, Phys. Rev. D, 37, 1237 
Sikivie, P. 1983, Physical Review Letters, 51, 1415 
Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, Astron. 

J., 131, 1163 
Taylor, G. B. 1996, ApJ, 470, 394 
Taylor, G. B. & Perley, R. A. 1993, ApJ, 416, 554 
Turok, N. 1996, Physical Review Letters, 76, 1015 
Vogt, C. & Enfilin, T. A. 2003, Astron. & Astrophys., 412, 373 
Vogt, C. & Enfilin, T. A. 2005, Astron. & Astrophys., 434, 67 
Weisskopf, M. C, Tananbaum, H. D., Van Speybroeck, L. P., & 

O'Dell, S. L. 2000, in SPIE Conference Series, Vol. 4012, SPIE 

Conference Series, ed. J. E. Truemper & B. Aschenbach, 2-16 
Wouters, D. & Brun, P. 2012, Phys.Rev., D86, 043005 
Wouters, D. & Brun, P. for the H.E.S.S. collaboration 2013, to 

appear in the proceedings of the Rencontres de Moriond, 

1304.0700 



